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Anomalous  surface  waves  including  Love  waves  and  phase  reversed  Rayleigh 
waves  have  been  observed  from  explosions  at  UTS  and  from  the  eastern  Kazakh 
areas  of  the  Soviet  Union.  In  this  report  we  use  a  linear  model  for  tectonic 
strain  release  to  estimate  the  amount  and  type  of  prestress  required  to  pro¬ 
duce  these  anomalies.  An  important  use  of  these  results  is  to  guide  the  input 
to  fully  nonlinear  simulations  of  an  explosion  detonated  in  a  prestressed 
environment. 

V  (Continued/ 


EOlTlON  or  »  NOV  •!  IS  OBSOLETE 


I  * 

. 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  RAGE  P»Ma  D* *• 


SS^i&^cs. 


f 

Unclassified 


UCUWTY  CLAMWWWIt  OF  This  gw 


A3STRACT:  (Continued) 

Tectonic  strain  release  adds  energy  to  the  explosion  component  of  the 
seismic  radiation  because  elastic  energy  storea  in  the  medium  is  released 
when  a  zone  of  weakened  material  strength  is  created  by  the  explosion. 

In  the  model,  this  zone  is  treated  as  if  it  were  a  spherical  cavity  with 
the  dimensions  of  the  explosion-produced  failure  zone.  This  effective 
cavity  radius  requires  Independent  constraints  which  can  be  provided  by 
nonlinear  simulations. 

General  conclusions  are  as  follows: 

In  the  presence  of  stress  concentrations,  tectonic  release  en¬ 
hances  high  frequency  radiation  in  preferred  directions. 

Compressive  stress  concentrations  reduce  body  wave  amplitudes 
while  tensile  stress  concentrations  amplify  body  waves. 

Long  period  tectonic  surface  waves  depend  only  on  the  average 
prestress  field  and  are  unaffected  by  stress  heterogeneity. 

The  tectonic  surface  waves  reduce  to  a  monopole  plus  a  quadru- 
pole  field  superimposed  on  the  explosion  monopole. 

Specific  conclusions  are  as  follows: 

Rayleigh  wave  reversals  or  factor  of  two  enhancements  (depend¬ 
ing  on  the  horizontal  prestress  being  compressive  or  tensile 
relative  to  the  hydrostatic  prestress)  can  be  obtained  with 
homogeneous  or  average  shear  stresses  of  about  50  bars 
assuming  that  the  effective  cavity  radius  is  equal  to  the 
explosion  elastic  radius. 

A  maximum  local  shear  stress  of  about  a  kilobar  is  required 
before  the  tectonic  component  of  body  waves  becomes  comparable 
in  size  to  the  explosion  body  waves.  The  average  prestress 
may  be  much  lower,  however. 

The  anomalous  body  wave,  Love  wave  and  Rayleigh  wave  radiation 
from  PILEDRI VER  can  be  simultaneously  explained  by  a  compressive 
stress  concentration  of  about  a  kilobar  to  the  northeast  of 
the  explosion. 
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SYSTEMS.  SCIENCE  AND  SOFTWARE 


!.  INTRODUCTION 


Rygg  (1979)  reports  the  observation  of  two  phase  reversed 
Rayleigh  waves  from  underground  nuclear  explosions  In  the  eastern 
Kazakh  area  of  the  Soviet  Union.  Patton  (1980)  reports  that  several 
additional  explosions  from  the  same  area  also  appear  to  have 
Rayleigh  wave  reversals.  Figure  1.1  shows  two  Rayleigh  waves  from 
eastern  Kazakh  (Shagan  River)  explosions  with  clearly  reversed 
polarities  recorded  at  Kabul,  Afghanistan.  Goforth,  Rafipour  and 
Herrin  (1982),  using  a  phase-matched  filter  analysis,  fird  phase 
reversals  and  time  delays  of  about  three  seconds  for  several  Shagan 
River  explosions.  Such  observations  have  been  reported  in  other 
areas  as  well.  Perhaps  the  best  documented  case  of  a  Rayleigh  wave 
reversal  is  from  the  NTS  exolosion  PILEDRIVER  (Toksoz  and  Kehrer, 
1972).  Rayleigh  wave  reversals  are  usually  accompanied  by  the 
generation  of  large  amplitude  Love  waves.  In  the  case  of 
PILEDRIVER,  a  body  wave  anomaly  also  occurred,  with  a  strong 
reduction  in  amplitude  to  the  Northeast  (Hadley  and  Hart,  1979). 
Since  body  wave  magnitudes  are  used  to  estimate  explosion  yields,  it 
is  very  important  to  understand  what  variations  in  body  wave 
amplitudes  may  accompany  surface  wave  anomalies. 

In  this  report  we  use  a  model  for  tectonic  release  which  can 
produce  Rayleigh  wave  reversals  and  large  Love  waves.  The  model 
assumes  that  when  ar,  underground  explosion  is  detonated,  a  shatter 
zone  (region  of  weakened  material)  is  formed  around  the  cavity.  If 
the  material  near  the  explosion  is  initially  in  a  strained  state, 
then  elastic  energy  stored  in  the  medium  is  released  in  the  form  of 
seismic  waves.  These  tectonic  surface  and  body  wayes  Interfere  with 
the  surface  and  body  waves  of  the  explosion,  and  will  cause  Rayleigh 
wave  reversals  if  the  tectonic  Rayleigh  wave  is  reversed  in  polarity 
and  exceeds  the  explosion  Rayleigh  wave  in  anp’Hude.  If  the 
horizontal  prestress  is  compressive  compared  to  the  hydrostatic 
prestress,  then  the  tectonic  Rayleigh  waves  will  be  reversed  in 
phase  compared  to  tne  monopole  explosion  Rayleigh  waves. 
Conversely,  a  horizontal  prestress  that  is  tensile  compared  to  the 
hydrostatic  prestress  causes  Rayleigh  wave  amplification. 
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Figure  1.1  Rayleigh  waves  from  two  explosions  in  the  eastern  Kazakh 

(Shagan  River)  area  of  the  Soviet  Union  which  exhibit  phase 
reversals.  The  t«o  events  were  recorded  at  station  KAAO  on 
18  August  1979  and  28  October  1979,  respectively. 
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The  explosions  whicn  snow  the  greatest  amount  of  tectonic 
release  are  those  in  granite  (Toksot  and  Kehrer,  1372).  For  our 
explosion  models,  we  therefore  use  two  spherically  symmetric  granite 
source  models  —  the  empirical  model  of  Mueller  ana  Murphy  (1971) 
and  a  numerical  model  for  partially  saturated  granite  (Bache, 
Barker,  Rimer  and  Cherry,  1980).  Re  add  the  wavefield  generated  by 
the  tectonic  release  model  to  these  monopole  explosion  sources  in 
order  to  determine  tne  effect  of  model  parameters  on  body  waves  and 
surface  waves. 

The  model  used  nere  to  simulate  tectonic  release  is  the  sudden 
creation  of  a  spnerical  cavity  in  a  prestressed  medium  (Stevens, 
1980).  This  provides  a  linear  approximation  to  the  complex 
nonlinear  processes  which  accompany  an  explosion  in  a  prestressed 
medium.  In  this  report  we  use  the  linear  model  to  estimate  the 
prestress  level  necessary  to  cause  Rayleigh  wave  reversals.  This 
estimate  may  be  used  to  provide  a  reference  level  for  fully 
nonlinear  simulations  of  an  explosion  in  a  prestressed  environment. 
Nonlinear  calculations  in  turn  may  be  used  to  provide  an  estimate  of 
tne  effective  “cavity"  radius  which  is  a  poorly  constrained 
parameter  in  the  linear  model. 

Day,  Rimer,  Cherry  and  Stevens  (1981)  performed  a  two-dimen¬ 
sional  non'inear  finite  difference  calculation  of  an  explosion  in  a 
prestressed  medium  where  the  prestress  consisted  of  overburden 
pressure  superimposed  on  a  75  bar  she'r  stress.  The  results  of  the 
calculation  snow  that  tne  effective  cavity  radius  (determined  by 
comparing  surface  wave  amplitudes  from  the  nonlinear  simulation  with 
surface  wave  amplitudes  from  the  linear  model  used  in  this  report) 
is  about  60  -  80  percent  of  the  elastic  radius.  Also,  the  waveforms 
of  long  period  surface  waves  found  in  the  nonlinear  simulation  are 
almost  identical  to  tne  waveforms  generated  using  the  linear  mooel. 

In  addition  to  finding  the  prestress  necessary  to  cause 
Rayleigh  wave  reversals,  we  also  estimate  the  amplitudes  of  the 
associated  Love  waves  and  the  variations  in  body  wave  amplitudes 
which  accompany  tne  reversals.  We  explicitly  include  the  effects  of 
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nearby  stress  concentrations  and  look  at  the  effects  of  the  stress 
innomogeneity  on  body  and  surface  waves. 

An  important  prediction  of  the  work  is  that  an  inhomogeneous 
prestress  field  should  have  a  much  stronger  effect  on  the  high 
frequencies  than  a  homogeneous  stress  field  would  have.  The  body 
waves  can  therefore  be  strongly  affected  by  inhomogeneous  stress 
fields,  especially  in  larger  explosions.  Also,  the  variation  in 
body  waves  can  be  highly  directional  with  a  strong  decrease  (or 
Increase)  of  body  waves  in  preferred  directions.  If  a  stress 
concentration  is  located  close  to  the  explosion,  the  body  waves  will 
be  most  strongly  affected  in  the  direction  of  the  stress 
concentration. 

Uniform  prestress  fields  in  general  should  have  a  much 
stronger  influence  on  explosion  surface  waves  than  on  the  body 
waves.  Assuming  that  the  region  with  minimal  shear  strength  extends 
to  the  elastic  radius,  Rayleigh  wave  reversal  may  be  obtained  with  a 
tectonic  prestress  of  about  SO  bars.  In  the  linear  model,  assuming 
that  shatter  zone  formation  is  instantaneous,  there  is  no  time  delay 
associated  with  the  reversal.  A  time  delay  would  exist  if  the 
shatter  zone  farmed  over  an  extended  period  of  time  (e.g.,  Stevens, 
1982) . 

Because  of  the  boundary  condition  of  vanishing  normal  and 
shear  tractions  at  the  free  surface,  the  only  truly  uniform 
prestress  fields  which  can  exist  are  linear  combinations  of  “strike 
slip"  and  uniaxial  prestress  fields.  The  tectonic  radiation  from 
explosions  in  these  prestress  configurations  has  negligible 
influence  on  steeply  descending  body  waves  from  the  explosion.  If 
the  prestress  field  is  not  quite  homogeneous,  then  the  prestress 
orientation  may  be  different  at  the  explosion  depth  than  at  the  free 
surface.  Some  other  prestress  orientations  can  have  stronger 
effects  on  the  explosion  body  waves.  The  results  of  the 
calculations  for  all  uniform  stress  fields  are  summarized  in  Table 
3.3  and  Figures  3.4  through  3.7. 
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Monuniform  prestress  fields  can  occur  in  many  different  ways. 
They  always  require  some  sort  of  dislocation  in  the  medium  since  a 
homogeneous  medium  without  dislocations  can  only  support  a  uniform 
stress  field.  One  way  in  which  an  inhomogeneous  stress  field  may 
form  is  by  faulting  under  the  influence  of  a  high  shear  stress 
field,  perhaps  when  the  material  is  deep  in  the  earth,  and  then 
removing  the  external  shear  stress  by,  for  example,  a  gradual  rise 
of  the  material  to  the  earth's  surface.  The  result  is  a  stress 
concentration  frozen  into  the  material  with  a  highly  inhomogeneous 
stress  field. 

Nonuniform  stress  fields  produce  more  complex  effects.  We 

analyze  the  body  and  surface  waves  when  explosions  occur  near  two 
types  of  stress  concentrations  —  a  point  static  dislocation  and  a 
center  of  compression  both  located  two  elastic  radii  from  the 
explosion  center.  The  results  are  similar  for  both  cases.  The 

stress  field  from  a  stress  concentration  varies  as  the  Inverse  cube 
of  the  distance  from  the  stress  concentration,  causing  a  strong 
variation  in  prestress  across  the  shatter  zone.  If  the  stress 

concentration  is  near  the  direction  of  the  body  wave  takeoff  angle, 
there  can  be  large  effects  on  the  explosion  body  waves.  Surface 

waves,  on  the  other  hand,  are  affected  only  by  the  average  stress 
field.  A  stress  concentration  which  causes  a  prestress  of  about  100 
oars  at  the  cavity  center  (800  bars  at  maximum  at  the  elastic 
radius)  can  generate  surface  waves  and  body  waves  comparable  to  the 
explosion  surface  and  body  waves.  Note  that  this  high  stress  field 
needs  to  exist  only  locally  and  at  depth  to  cause  body  wave 
anomalies.  A  summary  of  the  results  of  calculations  for  a  number  of 
com iguratlons  of  stress  concentrations  near  explosions  is  given  in 
Table  3. A  and  Figures  3.9  through  3.14. 

While  surface  wave  anomalies  have  been  widely  reported, 
reports  of  body  wave  anomalies  are  less  common.  However,  since  body 
waves  are  observed  from  only  a  small  part  of  the  focal  spnere  and 
may  be  strongly  affected  by  travel  path  effects,  body  wave  anomalies 
will  generally  be  difficult  to  detect.  A  clear  body  wave  anomaly 
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was  ccs.rved  from  the  explosion  PILCDPIVER,  however.  In  oritr  to 
estimate  whet  hinds  of  ,'nonal  ies  my y  be  exnectet  in  Soviet 

explosions,  ie  have  c-.r.ttrjcted  a  model  ol  the  prestress  a  id  which 

coaid  explain  the  PILEDRIVEP.  anomalies. 

The  sxplxsicn  P ILE0R1  Vt£R  ^ad  three  anomalous  fattures.  It, 
generated  large  Love  waves,  the  Rayleigh  waves  had  a  pronounced 
radiation  pattern  with  polarity  reversals  to  the  Northeast,  and  the 
uody  waves  were  gieatly  reduced  in  amplitude  at  stations  to  the 

Northeast.  We  can  explain  those  facts  simultaneously  by  assuming  a 
strong  coopressional  stress  concentration  to  the  Northeast  of  the 
explosion.  A  number  of  different  types  of  stress  fields  could 
explain  these  observations,  but  it  must  be  concentrated  tn  the 

Northeast  with  i  maximum  stress  of  about  a  xi lobar.  Rayleigh  wave 
reversals  are  caused  by  an  inward  motion  (parallel  to  the  surface) 
in  this  direction.  Small  body  nave  amplitudes  occur  near  a  node  in 
the  body  wave  radiation  pattern.  There  are  other  possible 
explanations  fo.’  the  PILEORiYER  cody  wave  anomalies,  Hadley  and 
Hart(1979)  for  examp’e,  suggested  mat  the  deep  structure  of  the 
climax  slock  was  rescansible.  The  very  large  surface  wave  anomalies 
observed  from  this  explosion,  however,  suggest  the  presence  of  a 
large  tectonic  stress  field  which  coulc  affect  the  body  waves  as 
well,  furthermore,  recent  stress  measurements  in  the  climax  granite 
(heuze,  e^al_. ,  i960)  found  the  nonhydrostetic  stresses  to  be  large 
and  highly  variable. 

Tectonic  strain  release  is  one  cf  the  few  siecnanisas  which  use 
explain  Rayleigh  wave  re/ersalt.  Masce  (19E1)  reviews  a  number  of 
hypotheses  whicn  have  oeen  made  fo-  this  end  ot.ier  surface  wave 

anomalies.  A  common  explanation  nas  aean  tnat  the  reversals  are 

caused  by  spall  slapdown  following  the  explosion.  However,  Day, 
Riser  and  Cherry  (1981),  show  that  if  momentum  is  to  he  conserved, 
then  spall  cannot  cause  Rayleigh  wave  reversals  and  that  spall  has 

only  minimal  effects  on  surface  wave  amplitudes  at  periods  longer 

than  ten  seconds.  Similar  problems  arise  with  other  passive 
mechanisms  for  producing  Rayleigh  wave  reversals. 


These  results  have  import!  c  consequences  for  yield  estimation 
ono  earthquake-explosion  discrimination.  It  is  possible  for 
Rayleigh  wave  reversals  to  occur  without  affecting  body  wave 
amplitudes;  however,  the  presence  of  surface  wave  anomalies 
increases  the  possibility  of  oody  wave  anomalies.  This  introduces 
®  a  cias  into  yield  estimates  which  may  be  difficult  to  remove.  Body 

waves  may  be  amplified  or  reduced,  or  both,  (at  different  receiving 
stations)  by  tectonic  release  although  this  effect  is  reduced  to 
some  extent  by  the  fact  that  the  explosion  body  wave  and  tectonic 
i  body  wave  are  not  exactly  in  phase,  and  so  will  not  cancel 

completely  at  high  frequencies.  Clearly  it  is  important  to  have  as 
broad  a  coverage  of  the  focal  sphere  as  possible.  In  the  case  of 
PILEORIVER,  strong  reductions  in  body  wave  amplitudes  occurred  over 
t  a  limited  range  of  azimuths,  but  our  models  show  only  a  small  effect 

on  body  wave  amplitudes  in  other  directions.  This  highly 

directional  behavior  is  characteristic  of  stress  concentrations. 

A  practical  method  for  removing  the  effects  of  tectonic 
release  would  be  to  apply  modern  invetse  theory  techniques  to  the 
body  and  surface  waves  simultaneously.  Cur  expansion  in  multipolar 
coefficients  lends  itself  well  to  this  type  of  analysis.  In  effect, 
this  procedure  would  be  an  automated  version  of  the  trial  and  error 
method  used  here  for  PILEDRIVER  whicn  would  consider  all  possible 
stress  fields  simultane  u'ly.  It  should  also  be  superior  to 
commonly  used  moment  l  Inversion  schemes  since  we  have  an 
estimate  of  the  frequency  dependence  of  the  composite  source  and  are 
not  restricted  to  a  monopole  plus  double  couple  source. 

The  effect  of  tectonic  release  on  discrimination  is  not  as 
severe  In  most  cases  as  the  effect  on  yield  estimates,  especially  if 
coverage  is  good.  Although  tectonic  release  adds  an 
“earthquake-line”  radiation  pattern  and  reduces  m^  (and  H$)  in 
certain  directions,  the  general  effect  is  to  add  energy  to  the 
explosion  body  and  surface  waves.  Since  the  explosion  and  tectonic 
body  waves  are  more  nearly  in  phase  at  longer  periods,  destructive 
t  interference  will  reduce  the  long  period  amplitude  relative  to  the 
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short  period  amplitude.  Furthermore,  Inhomogeneous  prestress  fields 
preferentially  add  high  frequency  energy  whicn  tends  to  increase  the 
corner  frequency.  Since  discrimination  methods  based  on  spectral 
analysis  such  as  the  Variable  Frequency  Magnitude  method  (e.g., 
Savino,  Archambeau  and  Masso,  1980)  compare  a  low  frequency 
amplitude  with  a  high  frequency  amplitude,  both  of  these  effects 
Improve  discrimination.  Nevertheless,  tectonic  release  can  cause 
discrimination  problems  in  some  cases.  Tensile  stress  fields,  for 
example,  may  increase  significantly  while  leaving 

unchanged,  thus  causing  a  more  “earthquake-like0  event. 

It  should  be  emphasized  that  in  this  report  we  have  considered 
only  tectonic  release  as  the  cause  of  anomalous  radiation  from 
explosions.  We  have  done  this  in  order  to  determine  the  nature  and 
size  of  prestress  fields  necessary  to  cause  the  observed  anomalies, 
not  to  rule  out  all  other  explanations.  It  is  possible,  for 
example,  that  the  Rayleigh  wave  reversals  from  PIlEDRIVER  were  due 
to  tectonic  release  in  a  uniform  stress  field  while  the  body  wave 
anomaly  was  caused  by  a  structural  feature.  So'  wa»e  anomalies,  in 
general,  do  require  much  higher  stresses  than  sonar'  wave  anomalies 
and  uniform  stress  fields  will  not  strongly  af-'ect  body  waves. 
Nevertheless,  we  find  nothing  inconsistent  with  tile  hypothesis  that 
stress  concentrations  may  exist  and  cause  variations  in  body  wave 
amplitudes  such  as  those  observed  from  PILEDRIVER.  This  must 
therefore  be  considered  as  a  possible  source  c'  bias  in  yield 
estimates. 


II.  A  MODEL  FOR  SEISMIC  RADIATION  FROM  AN  EXPLOSION 
IN  A  PRESTRESSED  *DIUM 


V  I 

i 

.•  >11 
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Almost  ail  nuclear  explosions  generate  some  "anomalous" 
radiation  in  the  fora  of  Love  waves  and  nonisotropic  Rayleigh  waves 
(e.g.  Toksoz  and  Kehrer,  1972).  In  some  cases,  however,  the 

anomalous  effects  can  De  very  large  and  actually  dominate  the 

obset .  t.iens.  There  are  a  number  o*  possible  explanations  for  this 

including  release  of  tectonic  prestress,  anisotropy  of  the  medium, 
nonsphericity  uf  the  sourue.  near  source  scattering,  slippage  on 
joints  and  faults,  and  induced  eartnquakes.  An  excellent  review  of 
the  suoject  can  ae  found  in  a  recent  paper  by  Masse  (1981).  All  of 
these  mechanisms  will  cause  the  generation  of  some  shear  waves  and 
cause  nonisotropic  radiation  patterns.  It  is  difficult,  however,  to 
explain  the  large  anomalies  of  FILEDR1VER  and  other  explosions  in 
which  the  Rayleigh  waves  actually  reverse  polarity.  Since  tha 

"anomalous"  source  wust  actually  excetd  the  direct  source  at  long 
periods,  it  seems  unlikely  to  result  from  a  passive  effect  such  as 
scattering  or  slippage  on  joints  and  faults.  Tectonic  release  can 
explain  Rayleigh  wave  reversals  since  it  adds  energy  to  the 
explosion  from  energy  stored  in  tne  e.rth. 

Hadley  and  Hart  (1979)  oointed  out  that  the  body  waves  from 
PILEOPIVER  were  anomalous  as  well  as  the  surface  waves.  In 
particular,  the  amplitudes  were  greatly  reduced  at  observation 
points  to  tne  Northeast.  Tectonic  release  from  stress 
concentrations  can  cause  this  kind  of  variation  in  the  body  waves. 

Vie  use  the  following  model  for  confuting  pndy  and  surface 
waves  from  an  explosion  in  a  prestressed  med'usn.  For  the  explosion, 
we  consider  two  models.  The  first  is  the  empirical  model  -f  Mueller 
and  Murphy  (1971)  In  which  the  explosion  is  described  ty  a  pressure 
pulse  of  the  fora: 

P  (t)  .  (P^  ♦  ?0(;)  H(t) 

applied  at  tne  elastic  radius  (the  point  bevend  which  motion  if 
elastic).  The  constants  ?0,  z,  PBC  are  determined  empirically 
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and  are  functions  of  yield,  material  properties  and  source  depth. 
For  PILEDRXVtFs  with  a  yield  of  60  kilotons  and  a  depth  of  470  meters 
in  granite,  these  constants  are  P0  »  130  bars,  ?oc  -  30  bars, 
a  «  15  sec-1.  The  elastic  radius  is  725  meters.  In  the  second 
source  model,  a  finite  difference  calculation  is  performed  for  the 
explosion  to  determine  the  reduced  velocity  potential  which  is  then 
used  to  synthesize  body  and  surface  waves  (Model  No.  469  from  Bache, 
8arker,  Rimer,  and  Cherry,  1380).  The  elastic  radius  for  the 
PIUDRIVER  finite  difference  calculation  is  870  meters. 

The  seismic  radiation  from  tectonic  release  is  modeled  using 
tne  solution  of  Stevens  (1980)  for  the  sudden  creation  of  a 
spherical  cavity  in  a  prestressed  medium.  This  solution  is  reviewed 
in  Appendix  1.  8y  “cavity"  we  mean  a  region  in  which  all  material 
strength  is  lost.  The  density  may  remain  constant.  A  problem  with 
this  model  is  how  to  determine  the  effective  cavity  radius. 
Tectonic  stress  will  be  released  thoughout  a  region  in  which  a 
substantial  amount  of  cracking  occurs.  This  region  will  be 
considerably  larger  than  the  explosion  cavity.  On  the  other  hand, 
*t  must  be  smaller  than  the  elastic  radius.  For  the  purposes  of 
this  study,  we  take  the  effective  “cavity  radius”  to  be  equal  to  the 
elastic  radius.  The  nuabers  we  find  will  therefore  be  the  lower 
limit  of  tne  actual  prestress  field.  A  smaller  cavity  radius  will 
require  a  larger  prestress  field  for  the  same  effect.  The  scaling 
relations  are  derived  in  Appendix  4.  It  is  approximately  true  that 
the  generation  of  seismic  waves  varies  as  the  product  of  the  stress 
field  and  the  cube  of  the  radius.  A  factor  of  two  decrease  in 
radius  requires  a  factor  of  eight  increase  in  stress  field  to 
produce  seismic  waves  of  the  same  amplitude.  This  is  not 
necessarily  true  for  body  waves,  however,  because  of  the  complexity 
of  tne  spectrum  at  short  periods. 

We  can  solve  for  the  seismic  radiation  from  a  prestressed 
snatter  zone  for  any  prestrass  field  (in  a  uniform  whole  space)  by 
performing  an  expansion  of  the  field  in  vector  spherical  harmonics. 
In  Aopendix  2  we  derive  the  expansions  for  all  uniform  stress  fields 
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and  for  two  types  of  stress  concentrations  —  a  point  center  of 
compression  and  a  point  (static)  dislocation. 

We  can  sake  synthetic  body  waves  and  surface  waves  from  these 
complex  sources  by  reducing  the  vector  solutions  to  scalar 
potentials.  Tills  conversion  is  given  In  Appendix  3.  The  potential 
solutions  can  be  expanded  in  cylindrical  potentials  and  Incorporated 
In  layered  space  codes  using  the  method  of  Bache  and  Harkrlder 
(1976)  for  body  waves,  and  the  method  of  Harkrlder  and  Archarabeau 
(1982)  and  Harkrlder  (1964)  for  surface  waves. 

Bache  (1976)  did  a  study  of  the  effects  of  tectonic  release  on 
P -waves  using  the  theory  of  Archambeau  (1972)  for  the  tectonic 
source.  He  concluded  that  tectonic  release  has  a  negligible  effect 
on  the  body  waves.  Archambeau' s  theory  was  valid  only  for  uniform 
stress  fields,  however.  Using  this  model,  Bache  concluded  that 
tectonic  release  could  not  have  a  significant  effect  on  body  waves 
without  generating  much  larger  surface  waves  than  are  observed  and 
without  assuming  much  higher  stresses  than  seemed  reasonable.  The 
situation  changes,  however,  when  Inhomogeneous  stress  fields  and 
stress  concentrations  are  considered.  The  following  analysis  will 
snow  tnat  concentrated  prestress  fields  can  affect  body  waves  more 
than  surface  waves.  It  Is  still  true  that  stresses  must  be  large, 
but  they  need  only  be  large  locally. 

As  shown  In  Appendix  1,  the  displacement  field  for  the 

tectonic  source  can  always  be  written  in  the  form: 

»  & 

U(X,  u)  *  Z  E  o2  !u)  N.n(«,x)  +  St  (u)  M^tu.x) 

Z« o  m*-i  “ 

*  Tta(u)  bn(u,i> 

The  coefflce.nts  o^,  sla,  depend  on  the  cavity  size  and 
type  of  prestress.  The  elgr.envectors  M.Q,  Sw  i,.m  represent 
toroidal,  spheroidal  shear,  and  spheroidal  compresslonal  waves 
respectively.  In  the  far  field  these  vectors  become: 


» 
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,1+1  e*V 


TT-^<i  +  1>  c^e.  t) 


Nln  *)  ”  'V'r8~  sim  (•»  *) 

hw  (••;)“  t lsrf  !»  <*•  e> 


where  the  vectors  Plm,  ^n,  ^  are  defined  in  Appendix  1  and 
depend  only  on  the  observation  coordinates  e,  t. 

The  far  field  displacement  spectrum  therefore  consists  of  sums 
of  terms  of  the  form  ana'll*,  «£_(“)/».  ■*£„(“)/“  multiplied 
by  angular  and  propagation  factors.  At  high  frequencies  the  spectra 
all  fall  off  as  w  .  At  low  frequencies,  the  function  o2n)(u)/u 
falls  off  as  a2-2.  The  functions  •»„(«)/».  t»_(<u)/b  are 

jwii  x,  m  _  y 

combinations  of  two  functions  which  fall  off  as  u  and  ■  "  ‘  at  low 
frequencies.  Thus  for  x  >  2  all  coefficients  vanish  in  the  limit 
b  >■  0.  At  sufficiently  low  frequencies,  any  source  must  therefore 
look  like  a  monopole  plus  a  quadrupole.  No  matter  how  complex  the 
source,  the  long  period  surface  waves  must  look  very  nearly  like  an 
explosion  plus  a  double  couple. 

For  uniform  stress  fields,  only  X-2  and  i- o  terms  exist.  For 
nonunifona  stress  fields,  however,  the  higher  order  terms  exist. 
These  terms  vanish  at  low  frequencies  and  peak  at  higher  frequency 
with  increasing  x.  The  frequency  of  the  peak  is  given  approximately 
by  FWJ(  •  xv/2tfi  for  SH  waves  and  by  (x+l)V/2ifi  and  (X-1)V/2tR  for 
the  twc  P-SV  waves  mentioned  before  where  V  Is  the  appropriate  wave 
velocity.  Inhomogeneous  stress  fields  can  therefore  cause  much 
larger  variations  in  the  body  waves  than  in  the  surface  waves, 
especially  for  large  explosions.  Higher  order  term  will  be 
important  when  V/2xR  <  fB  where  f  Is  the  dominant 
frequency  of  the  recording  Instrixsent.  For  PltEORIVER,  this  ratio 
is  about  1  Hx.  In  addition  to  adding  high  frequency  energy,  the 
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higher  order  multipoles  have  a  more  complicated  variation  with  angle 
so  inhomogeneous  stress  fields  can  give  highly  directional 
variations  in  body  waves. 

In  Appendix  2  we  derive  the  multipole  coefficients  for  two 
types  of  stress  concentrations.  In  the  following  section  we  will 
compute  synthetic  seismograms  and  body  and  surface  wave  radiation 
patterns  using  these  sources.  As  a  general  rule,  if  one  part  of  the 
“cavity"  is  more  highly  prestressed  than  another,  larger  amplitude 
tectonic  body  waves  will  be  generated  in  the  direction  of  the  higher 
prestress.  Depending  on  the  sign  of  the  prestress,  the  explosion 
body  waves  will  be  either  amplified  or  reduced  in  this  direction. 
Surface  waves,  on  the  other  hand,  will  be  sensitive  only  to  the 
average  stress  field  on  the  surface. 
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III.  THE  EFFECTS  OF  HOMOGENEOUS  AND  INHOMOGENEOUS  PRESTRESS 
FIELDS  ON  SURFACE  AND  BODY  WAVES 


In  this  section  we  mate  a  quantitative  comparison  of  the 
effects  of  different  types  of  prestress  fields  on  body  and  surface 
waves.  We  have  synthesized  seismograms  from  a  number  of  different 
prestress  fields  and  made  amplitude  comparisons  and  radiation 
patterns  for  each. 

Our  purpose  in  this  section  is  to  make  estimates  of  the  effect 
of  tectonic  release  and  the  magnitude  of  the  prestress  necessary  to 
cause  significant  changes  to  the  explosion  seismic  waves.  We  have 
used  earth  models  appropriate  for  PILEDRI VER  in  all  our  calcula¬ 
tions.  The  results,  however,  are  not  very  sensitive  to  the 
structure  and  should  be  applicable  to  most  regions  of  the  earth. 
The  structures  usea  are  given  in  Tables  3.1  and  3.2.  For  the  body 
waves,  the  source  region  is  a  PILEDRI VER  climax  stock  model  (after 
Bache  et.  al_.,  1975).  The  receiver  structure  is  a  standard  Eastern 
United  States  receiver  structure.  For  surface  waves,  we  used  the 
PILEORIVER  structure  over  the  Eastern  United  States  model  SI  of 
Bache,  Swanger  and  Shkoller  (198(1).  The  body  wave  calculation  used 
t*  ■  0.8.  The  body  wave  synthetics  were  made  using  a  WWSSN  short 
period  instrument  response.  Tne  upgoing  body  waves  were  reduced  by 
half  to  account  approximately  for  near  surface  attenuation  at  one 
hertz.  The  surface  wave  synthetics  were  made  through  an  LRSM  long 
period  instrument.  Ali  synthetics  were  made  using  a  cavity  radius 
equal  to  the  ftieller-Murphy  elastic  radius  of  725  meters.  This  is 
somewhat  smaller  than  the  finite  difference  elastic  radius  of  870 
meters.  The  receiver  point  was  taken  to  be  at  a  distance  of  4000 
kilometers.  This  requires  a  takeoff  angle  of  28  degrees  for  the 
body  waves.  The  Rayleigh  and  Love  wave  synthetics  include  only  the 
fundamental  mode. 

Synthetic  seismograms  for  the  explosion  source  are  shown  in 
Figure  3.1.  In  this  and  the  following  figures,  the  body  waves  are 
normalized  to  the  amplitude  of  the  Mueller-Murpny  explosion  body 
wave.  Surface  (Rayleigh  and  Love)  waves  are  normalized  to  the 
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TABLE  3.1 

CRUSTAL  HODELS  FOR  BODY  KAYE  CALCULATIONS 


Depth  (Km) 

Thickness 

a 

3 

P 

(Km) 

OCo/sec! 

(Km /sec) 

(qm/cn^) 

t 

SOURCE  REGION 

.05 

.05 

1.50 

1.00 

1.50 

.10 

.05 

4.50 

2.65 

2.55 

2.00 

1.90 

5.33 

3.15 

2.67 

• 

4.00 

2.00 

5.90 

3.25 

2.70 

CO 

CO 

6.00 

3.30 

2.80 

r. 

RECEIVER  REGION 

.6 

.6 

3.70 

2.16 

2.10 

C 

2.6 

2.0 

4.55 

2.54 

2.20 

4.0 

1.4 

5.60 

3.14 

2.65 

5.0 

1.0 

6.00 

3.30 

2.80 

G 


C 


« 
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CRUSTAL  MODEL  FOR  SURFACE  HAVE  CALCULATIONS 


Depth 

Thickness 

a 

6 

0 

0 

.05 

.05 

1.5 

1.0 

1.5 

10 

.10 

.05 

4.5 

2.65 

2.55 

20 

2.00 

1.90 

5.33 

3.15 

2.67 

200 

4.00 

2.00 

5.90 

3.25 

2.70 

300 

6.20 

2.20 

6.10 

3.30 

2.85 

400 

13.20 

7.00 

6.30 

3.40 

2.94 

1200 

19.00 

5.80 

6.40 

3.45 

3.00 

1500 

35.00 

16.00 

6.60 

3.60 

3.09 

2000 

eo 

30 

8.10 

4.50 

3.35 

2000 
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Mueller-Murphy  Body  Wave 


a  i  a  a  «  i  *  i  i  i  u  a  ii  i'i  «  a  »  j  ii  * 


Time  (sec) 


Finite  Difference  Body  Wave 


.2.3 

Time  (sec) 


Mueller-Murphy  Rayleigh  Wave 


» 


Explosion  Synthetic  Sei sinograms 


Figure  3.1.  Synthetic  seismograms  for  explosion  sources  without  tectonic 
release.  Top  is  Mueller-Murphy  empirical  body  wave  source 
synthetic.  Middle  Is  body  wave  from  Finite  Difference  Cal- 
t  culatlon  of  reduced  velocity  potential.  Bottom  is  synthetic 

Rayleigh  wave  from  the  Mueller-Murphy  source. 
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Mueller-fiurphy  Rayleigh  wave.  For  comparison,  the  finite  difference 
Body  waves  are  shown.  They  are  similar  In  shape  with  a  slightly 
larger  overshoot  at  the  end  because  of  the  higher  frequency  content 
in  the  finite  difference  source.  The  body  wave  amplitude  Is  greater 
by  a  factor  of  2.3.  The  surface  waves  are  Identical  In  shape,  but 
lower  in  amplitude  by  a  factor  of  0.8. 

3.1  HOMOGENEOUS  STRESS  FIELDS 

He  first  want  to  show  the  effects  of  tectonic  release  from 
explosions  In  spatially  uniform  stress  fields  with  different 
orientations.  Figure  3.2  shows  the  body  and  surface  wave  synthetic 
selsmograias  for  a  uniform  stress  field  of  one  bar  in  which  a 
spherical  cavity  Is  suddenly  formed.  The  radiation  pattern  from 
this  source  is  equivalent  to  a  strike  slip  double  couple  {with  a 
strike  of  -90').  The  body  waves  are  very  similar  to  the  boay  waves 
from  the  explosion.  The  main  difference  Is  a  slight  broadening  of 
the  waveform  (and  a  change  In  sign  at  this  azimuth).  The  surface 
waves  in  these  and  all  of  the  following  examples  except  for  the 
o13  stress  field  are  identical  except  for  changes  in  sign.  The 

sources  are  small  enough  that  the  surface  waves  seen  through  a  long 

period  instrument  are  Insensitive  to  source  properties  otner  than 
the  seismic  moment. 

Figures  3.4  to  3.7  sncw  the  body  and  surface  wave  radiation 
patterns  for  a  number  of  different  uniform  stress  fields.  The 
coordinates  used  are  shown  in  Figure  3.3.  The  results  are 

summarized  in  Table  3.3.  The  numbers  shown  are  the  amplitudes 
relative  to  the  ftieller-Murphy  P ILEDRIVER  source  for  a  unit 

prestress.  The  shaded  regions  are  the  regions  where  the  sign  of  the 
wave  is  reversed  relative  to  the  explosion.  For  surface  waves,  this 
is  a  well-defined  region.  For  body  waves,  It  Is  less  well-defined 
since  the  waveforms  can  vary  in  shape  as  well  as  in  sign,  especially 
near  nodes. 

The  oj2  (F'gwe  3.4)  stress  field  causes  a  four  lobed 
pattern  of  body  and  surface  waves.  This  stress  field  Is  much  more 
efficient  at  generating  surface  waves  than  body  waves.  (Body  waves 
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Figure  3.2.  Synthetic  seismograms  from  tectonic  release  in  a  c.,  st’ess 
field.  The  three  seismograms  are  the  P-wave,  Rayleigh  wave 
and  Love  wave  respectively.  The  shape  of  the  Rayleigh  and 
Love  waves  is  nearly  independent  of  the  source  and  is  very 
similar  for  all  of  the  following  examples. 
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Coordinate  systems  used  for  calculations  of  synthetic 
so'saograns  and  radiation  patterns.  The  Xj  axis  is 
out  of  the  paper.  The  Xg  axis  is  taken  to  be  $  *  0 
and  North. 
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TABLE  3.3 


> 
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BODY  AND  SURFACE  WAVE  AMPLITUDES  FOR 
EXPLOSIONS  AND  HOMOGENEOUS  STRESS  FIELDS 


<• 

Source  Type  P  { 

28*  takeoff) 

Rayleiqh 

Love 

EXPLOSION 

» 

Mueller-Murpny  (R«725ra) 

1.0 

1.0 

0.0 

Finite  Difference  (R=870ra) 

2.3 

.80 

0.0 

c 

UNIFORM  STRESS  FIELDS  (oars)  R*725.t, 

1 

0^  >  1  (strike  slip)  *(40 

1 

-42 

1 

17 

1 

3? 

’33  *  1  *•  ’ll  ■  ' ' 

1 

150 

1 

-32 

c 

(45*  Thrust) 

0p  «  -1  (Uniaxial) 

1 

'325 

l 

"42 

1 

35 

> 

op  «  1  (Vertical  Oip  Slip) 

1 

-106 

1 

1  800 

1 

1000 

’ll  *  -1  °22  c33  ’  2 

1 

-J7 

0 

(cylindrical) 
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RAYLEIGH  WAVE 


LOVE  WAVE 


Figure  3.4,  Radiation  patterns  For  bod/  and  surface  waves  for  a  ci? 

uniform  prestress  field.  Shaded  areas  indicate  jolarny 
opposite  tothe  explosion  in  this  and  following  figures. 
Amplitudes  are  relative  to  the  feller-murphy  explosion 
source.  Love  and  Rayleigh  wave  amplitudes  are  normalized 
to  the  explosion  Rayleigh  wave.  Body  waves  are  calculated 
at  one  hertz  and  have  a  23  degree  takeoff  angle.  Surface 
wave  amplitudes  are  computed  at  20  seconds. 
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3.5.  Radiation  patterns  for  uniform  prestress  field. 
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ire  observed  near  a  node  In  the  radiation  pattern).  The  Rayleigh 
waves  will  be  equal  to  the  explosion  Rayleigh  waves  with  a  prestress 
of  42  bars.  Thus,  Rayleigh  waves  will  be  phase  reversed  when  the 
prestress  exceeds  a 2  bars.  Love  waves  will  have  more  than  twice  the 
amplitude  of  the  Rayleigh  waves.  Body  waves,  on  the  other  hand, 
will  equal  the  explosion  body  waves  with  a  much  higher  prestress  of 
440  bars. 

Figure  3.5  shows  the  body  and  surface  wave  radiation  patterns 
for  a  stress  field  of  •  -  1,  -  1.  This  is  equivalent  to 
a  45  degree  thrust  fault  (radiation  pattern).  This  stress 
orientation  is  somewhat  more  efficient  at  generating  body  raves.  It 
still  requires  a  stress  of  150  bars,  however,  to  equal  the  explosion 
body  wave.  The  explosion  body  waves  are  anql if led  at  all  azlroutns. 
The  tectonic  Rayleigh  waves  are  very  large  for  this  orientation, 
requiring  a  stress  field  of  only  32  bars  to  equal  the  explosion 
Rayleigh  wave.  The  Rayleigh  waves  are  reversed  at  all  azimuths  for 
this  stress  field.  Love  waves,  again,  have  a  four-lobed  pattern  and 
have  about  the  same  amplitude  as  the  Rayleigh  waves. 

Figure  3.6  shows  the  radiation  pattern  for  uniaxial 
compression  (ojj  -  -1).  This  is  similar  to  the  45  degree  thrust 
orientation  except  for  the  lack  of  the  vertical  tensile  component. 
The  equivalent  source  Is  a  45  degree  thrust  plus  ar.  Isotropic 
compression.  This  stress  field  causes  a  reduction  in  the  body  waves 
and  a  reduction  (or  reversal)  of  the  Rayleigh  waves  at  all 
azimuths.  Again,  the  change  in  body  waves  is  fairly  small, 
requiring  325  bars  to  equal  the  explosion.  Rayleigh  wave  reversals 
require  a  prestress  about  40  bars  and  Love  waves  are  about  the  same 
amplitude  as  Rayleigh  waves.  It  would  be  difficult  to  distinguish 
these  two  cases  observational ly. 

Figure  3.7  shows  the  radiation  patterns  for  a  stress 
field.  This  is  equivalent  to  a  vertical  dip  slip  double  couple. 
This  orientation  generates  very  small  Love  and  Rayleigh  waves  with 
different  frequency  content  than  the  other  stress  orientations.  It 
is,  however,  a  fairly  good  generator  of  body  waves,  requiring  about 
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100  bars  to  equal  the  explosion  P  wave.  The  body  waves  have  two 
lobes  causing  amplification  In  half  the  directions  and  reduction  In 
the  rest.  This  type  of  stress  field  can  have  an  unpredictable 
effect  on  ratios.  Oependlng  on  the  direction  of  the 

observation  points,  it  could  cause  either  an  increase  or  a  decrease 
In  mjj  without  affecting  Hs-  The  body  wave  generation  of  the  dip 
slip  and  45  degree  thrust  prestress  orientations  Is  a  fairly  strong 
function  of  takeoff  angle  varying  as  sin  Z»  for  the  vertical  dip 
slip  and  cos  2»  for  the  45  degree  thrust.  The  amplitudes  become 
equal  at  «  »  22.5  degrees. 

Because  of  the  boundary  condition  of  vanishing  shear  and 
normal  tractions  at  the  free  surface,  uniform  and  stress 

fields  cannot  exist;  so  the  examples  for  "vertical  dip  slip"  and  "45 
degree  thrust"  orientations  are  not  physically  possible.  These  may 
exist,  however,  as  the  lowest  order  terms  of  an  Inhomogeneous  stress 
field.  The  radiation  from  tectonic  release  will  then  look  like  the 
radiation  from  one  of  these  "forbidden"  homogeneous  prestress  fields 
at  sufficiently  long  periods.  Also,  an  apparent  o33  stress  field 
appears  when  the  Isotropic  component  is  subtracted  out  from  a 
uniaxial  or  cylindrically  symmetric  stress  field. 

3.2  INHOMOGENEOUS  STRESS  FIELDS 

We  next  want  to  demonstrate  the  effects  of  stress  concentra¬ 
tions  on  body  and  surface  wave  radatlon  patterns.  The  results  are 
summarized  In  Table  3.4.  The  body  wave  radiation  patterns  are  for  a 
frequency  of  one  hertz  and  are  for  far  field  infinite  spa'e  (no 
reflections)  P-Waves.  The  surface  wave  rad1at1on  patterns  are  for 
20  second  Love  and  Rayleigh  waves. 

As  our  first  example  of  a  stress  concentration,  we  use  a 
fairly  simple  source  —  a  center  of  compression  located  one  radius 
from  the  surface  of  the  explosion  shatter  zone  (Figure  3.8).  The 
effect  of  this  stress  concentration  is  to  cause  a  highly  stressed 
region  on  a  small  part  of  the  cavity  surface.  For  the  example 
chosen,  the  stress  field  Is  eight  times  larger  at  the  surface  of  the 
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TABLE  3.4 


BOOY  AND  SURFACE  HAVE  AMPLITUDES 
FOR  INHOMOGENEOUS  STRESS  FIELDS 


* 


» 


C 


Angle  -  Stress 

Polarity 

p 

Rayleigh 

Love 

Concentration  to  Oown 

Center  of  Compression 

1 

11 

1 

Si 

1 

TO? 

30° 

- 

[150* ,  30*] 

none 

4- 

none 

all 

i 

1 

1 

104 

76 

46 

- 

[190*.  350*] 

[SO*.  120*] 

45° 

+ 

[80*,  100*] 

[-60* ,  60*] 

1 

1 

1 

152 

w 

IT 

[220*,  320*] 

[45*.  135*] 

60° 

4- 

[0\  180*] 

-45*.  45*] 

1 

1 

1 

163 

37 

27 

none 

[30*,  150*] 

90° 

4* 

ail 

[-30* ,  30*] 

Vertical  Dip  Slip 

1 

1 

1 

Dislocation 

unr 

IF 

IT 

• 

all 

[40°,  140'] 

45' 

+ 

none 

[-40°,  40'] 
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cavity  in  the  direction  of  the  stress  concentration  than  at  the 
center.  The  stress  field  falls  off  with  the  cube  of  the  distance 
from  the  center  of  compression.  All  of  the  examples  are  normalized 
to  a  stress  of  1  bar  at  the  cavity  center. 

The  radiation  patterns  are  quite  sensitive  to  the  location  of 
the  stress  concentration;  so  we  show  a  sequence  of  figures  with 
radiation  patterns  from  stress  concentrations  at  different  depths. 
Again,  the  shaded  areas  indicate  reversal  with  respect  to  the 
explosion  polarity.  For  body  waves,  as  mentioned  before,  this  is 
only  approximate  since  both  the  phase  and  amplitude  of  the  waves 
vary  with  angle.  If  the  stress  concentration  is  directly  under  the 
explosion,  then  the  radiation  is  azimuthally  symmetric.  The  body 
waves  are  all  decreased  in  amplitude,  the  Rayleigh  waves  are 
increased  in  amplitude,  and  no  Love  waves  are  generated. 

figure  3.9  shows  the  body  and  surface  wave  radiation  patterns 
for  a  stress  concentration  located  30  degrees  from  down  (at  an 
azimuth  of  minus  90  degrees  on  the  figure).  The  tectonic  body  waves 
are  reversed  in  polarity  relative  to  the  explosion  at  almost  all 
azimuths,  but  the  amplitude  is  much  greater  in  the  diiection  of  the 
stress  concentration  than  in  other  directions.  The  effect  of  this 
stress  field  when  added  to  the  explosion  is  to  selectively  reduce 
the  body  waves  in  this  direction.  It  is  also  a  good  generator  of 
body  waves  relative  to  surface  waves,  iilnety  bars  (at  the  center, 
7H0  maximum)  is  required  to  cancel  the  explosion  body  wave  although 
small  phase  differences  forbid  exact  cancellation.  Rayleigh  waves 
would  be  amplified  at  all  azimuths  by  this  stress  field.  It  is  a 
better  generator  of  Rayleigh  waves  than  of  Love  waves. 

Figure  3.10  shows  the  radiation  patterns  for  a  stress 
concentration  45  degrees  from  down.  Again,  the  body  waves  are 
substantially  reduced  in  the  direction  of  the  stress  concentration. 
The  Rayleigh  waves  are  amplified  in  most  directions;  however,  for  a 
narrow  range  of  azimuths  in  the  direction  of  the  stress 
concentration  (and  in  the  opposite  direction),  the  Rayleigh  waves 
are  reversed.  This  stress  field  also  generates  large  amplitude  Love 
waves. 
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Figure  3.9 


RadULion  patterns  for  tectonic  release  from  a  stress  con¬ 
centration  located  30  from  down. 
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Figure  3.10 


'•  Radiation  patterns  from  a  stress  concentration  at  45*. 
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Figure  3.11  shows  the  radiation  pattern  for  a  stress 

•  concentration  located  60  degrees  from  down.  This  configuration 
still  causes  a  body  wave  reduction  in  the  direction  of  the  stress 
concentration.  It  is  less  efficient  than  the  previous  examples, 
however,  and  directions  away  from  the  stress  concentration  see  an 

t  amplification  of  the  body  waves.  There  is  a  large  region  in  which 

Rayleigh  waves  are  reversed.  At  some  azimuths,  they  are  amplified. 

Figure  3.12  shows  the  radiation  patterns  for  a  stress 

concentration  at  90  degrees  at  the  same  depth  as  the  center  of  the 
cavity.  Now  body  waves  are  amplified  at  all  azimuths.  Rayleigh 

waves  are  reversed  at  all  but  a  narrow  range  of  azimuths.  This 
stress  field  is  a  relatively  efficient  generator  of  Love  and 

Rayleigh  waves. 

*  Figure  3.13  shows  synthetic  seismograms  at  three  observation 

points  —  azimuths  of  minus  90,  0,  and  90  degrees  with  a  center  of 
compression  at  an  azimuth  of  minus  90  degrees  and  *5  degrees  from 
down  (radiation  pattern  of  Figure  3.10).  At  minus  90  degrees,  the 

i  tectonic  body  wave  has  almost  the  same  shape  as  the  explosion  and 

will  therefore  tend  to  reduce  the  body  wave.  At  90  degrees,  the 
wave  is  amplified.  At  0  degrees,  the  waveform  is  more  complex  and 
does  not  clearly  increase  or  decrease  the  explosion  body  wave. 

(  As  a  second  example  of  a  stress  concentration,  we  use  a  point 

static  dislocation  located  two  radii  from  the  cavity  center.  The 
dislocation  is  located  at  an  azimuth  of  minus  90  degrees  as  with  tn* 
center  of  compression  examples.  It  has  a  vertical  dip  slip 

C  orientation  and  is  located  45  degrees  from  down.  The  explosion  then 

is  in  a  compressive  quadrant  of  the  double  couple  radiation 

pattern.  Again,  the  stress  field  falls  off  with  the  cube  of  the 
distance.  The  radiation  patterns  for  this  source  are  shown  in 

C  Figure  3.14.  The  radiation  resembles  that  of  the  center  of 

compression.  It  is  slightly  more  efficient  at  generating  surface 
waves  and  slightly  less  efficient  at  generating  body  waves  than  the 
center  of  compression.  The  dislocation  at  45  degrees  has  radiation 

f  patterns  and  amplitudes  much  like  the  center  of  compression  at  60 

degrees. 
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Figure  3.11 


.  Radiation  patterns  from  a  stress,  concentration  at  E0\ 
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TIME  (sec) 

Synthetic  Body  Wives  From  Center  of  Compression  at  45' 

Figure  3.13.  Synthetic  body  wave  seismograms  at. different  azimuths.  The 
stress  concentration  Is  located  45 ‘from  down  at  $  *  *  90‘. 
The  wave  forms  cnangs  in  shape  with  azimuth  and  reduce  the 
body  wave  in  the  direction  of  the  stress  concentration. 
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Figure  3.14.  Radiation  patterns  from  a  vertical  dip  slip  static 
dislocation  at  45°. 
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In  the  preceding  example,  notice  that  the  stress  field  was 
»  caused  by  a  vertical  dip  slip  dislocation,  but  the  radiation 

generated  by  tectonic  release  In  this  stress  field  Is  very  different 
from  that  found  previously  for  a  “vertical  dip  slip*  homogeneous 
prestress  field.  It  Is,  in  fact,  very  efficient  at  generating 
I  surface  waves  and  has  much  the  same  effect  as  the  stress  field  of 

the  point  center  of  compression. 

It  Is  clsar  from  these  examples  that  the  relative  effects  of 
'tress  concentrations  on  body  and  surface  waves  Is  highly  dependent 

i.  on  the  location  of  the  stress  concentration.  The  body  waves  csn 

have  complex  radiation  pattens  —  patterns  which  could  not  occur 
with  a  simple  double  oouple  source,  it  is  approximately  true  that 
the  generation  of  body  waves  is  a  function  of  the  distance  from  the 

S  stress  concentration  to  the  surface  of  the  cavity  (shatter  zone) 

while  the  generation  of  surface  waves  is  a  furctlon  of  the  the 
distance  from  tne  stress  concentration  to  the  center  of  the  cavity. 
The  effect  on  body  waves  will  therefore  be  stronger  compared  to  the 

j.  effect  on  surface  waves  if  the  stress  concentration  is  moved  closer 

to  the  explosion.  This  is  especially  true  for  large  explosions.  As 

pointed  out  in  the  last  section,  the  higher  order  multipoles  add  In 
seismic  energy  at  and  above  the  comer  frequency.  If  the  comer 
frequency  is  much  higher  than  the  observation  frequency,  these 
effects  will  not  be  observed.  At  frequencies  higher  than  the  comer 
frequency,  these  effects  will  be  very  lnoortant.  The  effect  on 
surface  waves  will  usually  be  Indistinguishable  from  a  quadrupole 
plus  isotropic  source,  tody  waves  will  generally  be  Increased  In 
directions  corresponding  to  tensile  stresses  on  the  sphere  and 
decreased  in  directions  of  compressive  stress  on  the  sphere.  The 
stresses  necessary  are  higher  than  the  stresses  for  uniform  stress 
fields  for  body  waves  o*  the  same  amplitude.  For  example,  with  a 

c  center  of  compression  at  45  degrees  from  down.  It  requires  a  auxlnum 

stress  of  300  bars  to  equal  the  explosion  body  wave  conpare'  to  ZOO 
bars  for  a  “45  degree  thrust"  stress  orientation.  This  region  of 
high  stress  Is  quite  small,  however,  and  decreases  by  a  factor  of  27 
O  across  the  cavity. 
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IV.  A  TECTONIC  SOURCE  MODEL  FOR  PUECRIVER 


The  explosion  PILEDRI VER  emitted  highly  anisotropic  surface 
waves  and  body  waves.  Toksoz  and  Kenrer  (1972)  were  able  to  natch 
the  Rayleigh  wave  radiation  patterns  with  the  superposition  of  a 
strike-slip  double  couple.  In  their  moael,  the  amplitude  of  the 
Rayleigh  waves  created  by  the  double  couple  exceeded  those  of  the 
explosion  by  a  factor  of  3.2.  This  number  can  be  reduced  by  using  a 
different  orientation  of  the  double  couple  such  as  a  45  degree 
thrust  orientation,  but  it  is  still  a  very  large  effect.  The 
observed  Rayleigh  wave  polarity  was  reversed  at  azimuths  to  the 
Northeast. 

Hadley  and  Hart  (1979)  compared  the  body  waves  from  PILEDRIVER 
with  the  body  waves  from  one  explosion  JORUM  and  found  a  substantial 
variation  in  the  amplitude  of  the  body  waves  with  station  azimuth. 
In  particular,  the  body  waves  in  the  Northeasterly  directions  were 
reduced  by  as  much  as  a  factor  of  six.  Table  4.1  shows  the  body 
wave  amplitude  data  for  JORUM  and  PILDERIvER  and  the  ratio  of  the 
amplitudes  scaled  with  the  cube  root  of  the  yields.  The  data  is 
from  two  Sierra  Seophysics  reports  (Hadley,  1979,  and  Hadiey  and 
Hart,  1979). 

The  fact  that  the  Rayleigh  wave  reversals  are  in  the  same 
direction  as  the  body  wave  anomaly  suggests  that  they  both  have  the 

same  cause.  The  purpose  of  this  analysis  is  to  assume  that  the 

anomalies  are  due  to  tectonic  release  and  to  see  what  magnitudes  and 

orientations  of  stress  fields  are  necessary  to  match  the 
observations.  The  solution  is  not  unique,  but  there  are  some  strong 
constraints  on  the  problem.  Since  the  body  wave  variation  is  large, 
the  prestress  field  cannot  be  too  efficient  at  creating  surface 
waves  relative  to  body  waves.  A  “strike-slip*  (oj-)  prestress, 
for  example,  would  generate  very  large  surface  waves  without 
substantially  affecting  the  body  waves.  A  second  constraint  is  that 
the  Rayleigh  waves  are  reversed  and  the  body  waves  are  reduced  in 
the  same  direction.  This  means  that  the  prestress  must  be 
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TABLE  4.1 

BODY  WAVE  AMPLITUDES  FOR  EXPLOSIONS  PILEDRIYER  AND  JORUM 


Station 

Distance 

(decrees ) 

Azimuth 

(deqrees) 

J-Ajep 

(mu) 

P-Amp 

(mu) 

Scaled  PILEDRIYER 

JORUM 

SCP 

30.0 

71.2 

823 

69 

.214 

000 

32.4 

70.0 

745 

46 

.158 

COL 

33.2 

336.1 

1074 

343 

.816 

WES 

34.7 

67.0 

463 

34 

.186 

RIP 

39.9 

258.7 

1097 

548 

1.28 

KTG 

56.8 

23.7 

594 

91 

.391 

AXU 

59.9 

28.3 

365 

91 

.637 

NNA 

61.6 

134.8 

422 

171 

1.04 

ARE 

58.2 

133.1 

1268 

446 

.898 

tol 

81.3 

46.0 

685 

91 

.339 

STU 

31.7 

32.9 

251 

69 

.702 

SHK 

83.9 

309.1 

868 

114 

.335. 

NAT 

86.1 

99.8 

548 

183 

.853 
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compressive  over  a  large  part  of  the  focal  sphere  In  the 
*  Northeasterly  direction.  A  “45  degree  thrust*  stress  field,  for 

example,  causes  a  reduction  In  the  surface  waves  but  an 

amplification  of  the  body  waves.  In  fact,  the  only  “fault-1 1ke“ 
uniform  stress  field  which  could  satisfy  these  conditions  Is  a 
t  “strike-slip*  stress  field  oriented  approximately  North-South  with  a 

compressive  lobe  to  the  Northeast  and  then  rotated  down  to  the 
Northeast  so  that  there  Is  a  larger  effect  on  tne  body  waves. 


There  Is  no  reason  to  restrict  ourselves  to  uniform  stress 
fields.  In  fact,  the  stress  field  almost  certainly  Increases  with 
depth.  The  major  observational  restriction  Is  that  the  stress  field 
be  large  and  compressive  to  the  Northeast.  He  will  therefore  try  to 
match  the  body  and  surface  waves  with  a  stress  concentration  In  the 
•  lower  Northeast  quadrant.  From  the  figures  3.9  to  3.14  the  location 

of  the  compresslonal  stress  concentration  can  be  determined.  If  the 
stress  concentration  is  located  too  low  (0-45  degrees  from  down) 
then  there  is  a  large  effect  on  the  body  waves,  but  an  actual 
amplification  of  the  surface  waves  and  no  significant  Rayleigh  wave 
reversal.  If  the  stress  concentration  Is  at  the  same  depth  as  the 
center  of  tne  cavity  (90  degrees  from  down),  then  the  Rayleigh  waves 
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can  be  reversed,  but  there  is  no  reduction  In  the  body  waves.  At  60 
degrees,  nowever,  the  effect  is  right.  The  negative  Rayleigh  wave 
lobe  is  large  enough  to  cause  a  reversal  when  superimposed  on  the 
explosion  source.  Body  waves  are  reducid  preferentially  in  the 
Northeasterly  direction.  The  best  linear  combination  with  the 
Hieiler-Hurphv  explosion  source  is  a  center  of  compression  at  60 
degrees  from  down  at  an  azimuth  of  35  degrees  from  North  with  a  200 
bar  stress  field  at  the  cavity  center.  This  leads  to  a  maximum 
prestress  of  1.6  kilobars  on  the  “cavity*  surface. 

The  radiation  patterns  for  this  source  are  shown  in  Figure 
4.1.  The  Rayleigh  wave  radiation  pattern  agrees  very  well  with  the 
data  given  by  Toksoz  and  Kehrer  (1972,  p.150).  In  fact,  the  match 
Is  better  than  their  fit  with  a  superimposed  strike-slip  double 
couple. 
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Figure  4.1.  Radiation  patterns  for  composite  PILEDRIVER  source.  Marks 
«"  ^  «ve  f19ure  are  PILEDRIVER  observations  divided  by 
JORUM  observations. 
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Rivers  and  Von  Seggern  (1981)  did  a  moment  tensor  Inversion 
for  the  PILEDRIV8R  source  using  Rayleigh  waves  and  love  waves 
without  restricting  the  "double  couple"  direction  to  be 

strike-slip.  The  Love  wave  radiation  patterns  for  the  resulting 
source  are  very  similar  to  the  Love  wave  patterns  in  Figure  4.1. 

The  body  waves  shown  are  for  a  frequency  of  1  Hertz  and  a 
takeoff  angle  of  28  degrees.  The  data  points  shown  are  the 

amplitudes  .f  P ILEDRIVER  divided  by  the  amplitudes  of  JORliM  at  the 
stations  listed  in  Table  4.1.  The  takeoff  angles  actually  vary 
between  IS  degrees  and  30  degrees;  so  the  stations  farther  away 
should  have  a  more  uniform  radiation  pattern  than  Is  shown  on  the 
figure.  Agreement  is  quite  good.  The  one  bad  datum  is  the  low 
reading  at  station  SHK  (at  azimuth  of  310  degrees). 

Figure  4.2  shows  synthetic  seismograms  from  the  composite 
source  along  with  their  relative  amplitudes.  The  conditions  are  the 
same  as  for  the  synthetics  in  the  last  section  except  that  we  used  a 
value  for  t*  of  1.05.  This  value  was  obtained  by  Bache,  et  aJU 
(1975)  for  the  path  from  NTS  to  Alaska.  The  observed  amplitudes  are 
close  to  the  syntnetlcs.  The  Northeasterly  waveforms  are  more 
complex  aue  to  interference  of  the  explosion  and  tectonic  body 
waves.  The  synthetic  amplitude  at  0G0  Is  too  large  Indicating  more 
complete  destructive  Interference  than  Is  in  our  model. 

Other  types  of  stress  concentrations  to  the  Northeast  give 
similar  results.  From  the  calculations  In  the  iast  section,  a 
vertical  dip  slip  dislocation  located  45  degress  from  down  has 
almost  the  same  effect  as  the  center  of  compression  at  60  degrees. 
Again,  with  this  source  we  can  approximately  match  the  surface  and 
body  waves.  The  maximum  stress  on  the  surface  Is  about  1.5  kllobars. 

As  mentioned  before,  the  maximum  prestress  on  the  cavity 
surface  from  point  stress  concentration  is  about  1500  bars.  In 
fact,  the  stress  field  does  not  need  to  be  this  large.  A  somewhat 
more  extended  stress  field  can  also  cause  the  observed  effects. 
Figure  4.3  shows  the  radiation  patterns  with  a  stress  concentration 
at  45  degrees  from  down  (maximum  stress  on  cavity  800  bars)  at  an 
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Figure  4.2.  Body  wave  synthetic  seismograms  made  using  PIIEDRIVER 

composite  source.  The  complexity  of  the  source  causes  a 
variation  in  waveform  shape  with  takeoff  angle  and  azimuth. 
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azimuth  of  40  degrees  and  a  uniaxial  prestress  of  100  bars  along  the 
same  azimuth.  Again,  the  agreement  is  good,  but  the  maximum  stress 
field  Is  only  about  SSO  bars. 

The  crucial  condition  in  this  analysis  is  that  the  lower 
Northeast  quadrant  be  under  substantial  compressive  prestress.  The 
surface  waves  are  very  Insensitive  to  the  type  of  prestress  and  are 
really  only  affected  by  the.  monopole  and  quadrupole  components.  The 
body  waves  are  more  sensitive  to  details  of  the  stress  field.  From 
these  examples  It  Is  obvious  that  we  cannot  determine  the  exact 
magnitude  of  the  stress  field.  To  within  about  a  factor  of  two, 
however,  the  maximum  stress  should  be  approximately  a  kilobar. 
Furthermore,  the  stress  needs  to  be  high  only  locally  and  at  some 
depth. 

In  this  section  we  nave  demonstrated  that  the  general  features 
of  the  explosion  PILEDRIVER  —  the  Rayleigh  wave  reversals,  large 
amplitude  Love  waves,  and  asymmetric  body  wave  radiation  pattern  — 
can  be  explained  If  there  Is  a  compressive  stress  concentration  of 
at  least  a  kilobar  to  the  Northeast  of  the  explosion.  We  have  not 
attempted  to  make  a  detailed  model  of  the  stress  field  or  to  explain 
small  details  of  the  seismograms.  In  principle.  It  should  be 
possible  to  perform  an  Inverse  calculation  to  determine  the 
multipole  coefficients  and  the  Initial  prestress  field.  This 
inversion  could  be  done  using  body  waves  and  surface  waves 
simultaneously.  The  problem  Is,  of  course,  under-determined  since 
we  have  only  limited  coverage  of  the  focal  sphere,  but  would  provide 
a  method  to  describe  the  class  of  models  which  fit  the  data  In  a 
single  calculation  as  opposed  to  using  a  trial  .and  error  method  of 
forward  modeling. 


46 


SrxrtMt,  JC/CNCC  AN O  AOATNAAK 


) 


» 


» 


* 


* 


* 


t 


I) 


i 

■  !• 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 

Rayleigh  wave  reversals  nay  be  obtained  if  an  explosion  occurs 
in  a  medium  with  a  prestress  greater  than  about  50  bars,  with  the 
exact  amount  depending  on  the  orientation  of  the  stress  field, 
assuming  that  the  effective  shatter  zone  radius  Is  approximately 
equal  to  the  elastic  radius.  The  size  of  the  effective  radius  Is 
critical  since  the  long  period  surface  wave  amplitude  is  propor¬ 
tional  to  the  cube  of  the  -adius.  The  best  way  to  determine  the 
effective  radius  Is  by  performing  nonlinear  calculations  of 
explosions  in  prestressed  media  and  comparing  it  with  the  linear 
result  (Day,  et  ah ,  1081). 

An  Important  fact  in  light  of  recent  observations  is  tn't  the 
linear  model  used  here  does  not  produce  time  delays.  Any  time 
delays  must  result  from  tectonic  release  over  an  extended  period  of 
time.  The  only  way  to  determine  if  time  delays  arc  present  and  the 
extent  of  these  delays  without  making  ad  hoc  assumptions  is  by 
performing  full  nonlinear  calculations  of  explosions  with  tectonic 
release  using  realistic  material  models. 

Surface  waves  are  very  insensitive  to  stress  heterogeneity. 
Only  the  average  stress  field  is  important  for  long  period  Rayleigh 
wave  excitation  although  the  orientation  of  the  average  stress  field 
s  very  important. 

Explosions  near  stress  concentrations  may  have  significantly 
altered  body  waves.  These  effects  can  be  strong  in  some  cases  and, 
in  general,  will  be  highly  directional.  Local  stresses  of  about  a 
kl lobar  are  required  to  Influence  the  body  waves;  however,  these 
stresses  need  be  this  large  only  locally  and  at  depths  below  the 
explosion.  The  average  stress  may  be  much  lower  than  the  local 
stress  fields. 

We  have  used  the  model  for  tectonic  release  to  model  the 
explosion  PtlEDRIVER.  If  the  body  and  surface  wave  anomalies  are 
due  to  tectonic  release,  then  the  observations  require  that  the 
average  stress  field  nust  be  compressive  1r  the  northeast-southwest 
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direct? dr.  The  bouy  wave  anomaly  and  surface  wave  anomaly  can  be 
explained  simultaneously  py  a  stress  concentration  below  and  to  the 
northeast  of  the  explosion. 

Tectonic  release  causes  problems  for  discrimination  and  yield 
estimation  by  causing  variations  in  surface  and  body  wave  amplitudes 
at  different  locations.  The  most  important  step  (n  removing  these 
problems  is  to  have  as  wide  a  coverage  of  stations  as  oossible. 
Moment  tensor  inversion  is  a  good  way  to  remove  the  effects  'of 
tectonic  release;  however,  a  cylinarically  syimetric  stress  field 
canr.ot  be  removed  this  way  using  surface  waves  alone. 
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APPENDIX  I 

SOLUTION  FOR  THE  CREATION  OF  A  SPHERICAL 
CAVITY  IN  A  PRESTRESSED  MEDIUM 
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SOLUTION  FOR  THE  CREATION  OF  A  SPHERICAL 
CAVITY  IN  A  PRESTSESSEO  MEDIUM 

An  Integral  equation  for  the  displacement  field  caused  by 
tractions  T  (u_)-n  applied  to  the  inside  of  a  cavity  of  any  shape  can 
be  written  (in  the  frequency  domain)  as: 


(x,  »)  .  y"u  .  T  (G)  •  ndA  -  j  G  .  T  (u) 


ndA 


(1) 


where  G  (x,  «)is  any  frequency  domain  elastic  Green's  tensor, 
and  Z  is  the  cavity  surface.  If  the  cavity  is  created  suddenly  in  a 
prestressed  medium,  the  apparent  tractions  are 

T  (u)  •  n  .  -  i.  T  (u»)  •  J5 
where  T  (u*)  is  the  initial  prestress. 


For  a  spherical  cavity  equation  (I)  is  separable.  He  define 
the  complete  set  of  vector  spherical  harmonics  (after  Ben-Aienahem 
and  Singh,  1S68): 

£ta  *  *r  YW®’ 


Au*l)  ^  *e  +  %  sine  7? 

Ytm<®*  « 

c 

r  * 

£ta  “  [  e»  sTne^f-  %  75  J  1 

tm  <*•  *> 

(2) 

* 

Ytm  <*•  *>  ’  pln  <C0S8> 
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Then  If  we  expend  the  Green's  tensor  end  Its  derlvetives,  the 
displacement  and  the  prestress  In  terms  of  these  functions,  the 
orthogonality  relations: 


) 


I 


t 


» 


» 


z 


c 


* 


» 


/  V  ? 


and 


i'm' 


,  dSl 


dn*  /?£m  •  dn-/!tm*  £ J'm* 

*  / Ilm  *  dfl*/V  £*■•  dn' 


d  15 «  0 

(3) 


^mm1 


4t  i  1+  m  1 

•  rrrr  i  rrsfr 


used  at  the  cavity  surface  (X  -  X  )  reduce  equation  (1)  to  a  set 
of  algebraic  equations.  Oetalls  of  the  solution  may  be  found  In 
Stevens  (1980).  Given  the  tractions  In  the  form: 


T(u*)  •  n  .  £  £  0,  J,  P5o{»,  t )  ♦  Oj  Ml  *  1)  B*n(o,  t) 


" .  im  lin 
lm o  HW-i 


(4) 


the  displacement  may  be  expressed  in  terms  of  the  elastic 
eigenvectors: 


•O  i 

u  (x,  .)  -  £  £  °il5  («)  %'*•  •)  *  6tm^“) 

~  ~  i»0  ow-i 

*  W-)  (5) 

The  coefficients  ala,  s?a,  T  determine  the  displacement 
field  at  any  point  outside  the  cavity.  The  vector  eigenfunctions 


1 

i 
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Him*  4m*  km  rePr«ent  toroidal,  spheroidal  shear,  and 
spheroidal  congressional  waves  respectively.  They  are  defined  by: 

%-yy)  '/iUM!  ctm(9.  t) 

\  (y) 

%  *  iU*  k  -7- 

(y)  , - 

♦(f£(y)*— j — )  8^(4,  ()  (6) 

f,  (x)  - 

km  "  {i  <x>  %<«•  *)  *  J  AuTl)  BZn(4,  *) 

where  x  .  K^r,  y  .  Ksr,  and  fz  .  h'?>,  the  spherical  Hankel 
function  of  the  second  kind. 
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APPENDIX  2 

VECTOR  MUITIPOLE  COEFFICENTS  FOR  PRESTRESS  FIELDS 


I 


» 


I 


I 


% 


I 


II 


In  order  to  use  the  preceedlng  solution  It  Is  necessary  to 
nave  the  Initial  stress  field  expanded  in  vector  spherical  harmonics: 


T  (u«) 


E 

l  .0 


l 

E 


Dim  r_m 


♦  D^yiTiTTr  a) 

The  coefficients  D1^,  0^,  D^  may  always  be  found  by  using 
the  orthogonality  relations  of  the  vector  harmonics,  but  In  some 
cases  they  may  be  found  analytically.  A  uniform  stress  field 
(cijj  »  constant)  may  be  written  in  the  form: 


T(u’>  *  "  -  A00  loo  *  J2  *2m  (%+,-r  v 


(2) 


with 

Aoo  "  T  (ffll  *  °22  *  ’33* 
and  defining  the  deviatorlc  stress  tensor 

®ij  "  a1j  '  {1j  Aoo 
The  remaining  coefficients  are: 

*22  *  “ir  (®U  -  ’22)  *  5T  °12 


•1 
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A21  *  T  °13  +  4r 


The  coefficients  in  equation  (1)  are  then 


DSo  “  Aoo 


°L  ‘  \  A2m 


All  other  coefficients  vanish. 

The  negative  m  values  may  always  be  determined  from  the 
positive  m  values  by  using  the  fact  that  the  field  is  real.  This 
requires  that: 


(ff^a  is  the  complex  conjugate  of  D^) 


Two  other  prestress  fields  which  may  be  expanded  analytically 
are  a  point  center  of  compression  or  a  point  dislocation  located 
along  the  Z  axis  (see  Figure  AZ-1). 

The  stress  field  for  the  center  of  compression  may  be  easily 
found  by  using  the  fact  that  the  (static)  displacement  field  can  be 
written 


"  *  -  J?  «r  ■  c  ?  F 
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Compression 


I 


Figure  A2.1.  Coordinates  used  for  computation  of  stress  fields 
from  center  of  compression  and  static  dislocation. 
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Then  expanding  by: 


■  £  -T+r  P,  (cos#) 
1*0  L1  1 


;  -  S  e  -k  V C0S9) 

~  ImO  L 


The  tractions  due  to  any  displacement  field  of  the  form: 


!  *•  «  -  £  W*>  !tm  *  Wr)  !w  +  f3tm  St™  (9) 


are  given  by 


x(u>.er -  g  -yro^T)  %*) 


+  2“  -JT2  ptm  (•.  *) 


u  ^  -  ~S  j 

/,f3tm  f3im\  -  ,a 

BpiF - —j  Etm(#’ 


!w <»• « 
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tensor  (Ben-Menahen  and  Singh*  1968)  and  applying  the  stress 
operator  to  equation  (14),  the  stress  field  can  be  written: 


I  (u*)  '  ®r 


l  »C  lita-Z 


■  T  Vj,)  •  er  (fi  •  T  (N^)  •  l)/U  *  1H2Z  *  3)  1 

'  '  /1C) 


(15) 


where  the  external  eigenvectors  are: 

r*  l~l  U(  l  *  ij  8..  -  (!♦  1)  P 


Nta  -  r-  JvM  i  -  U  ha'  t;  ran 


c-  r-Zl 
F .  «  r  — 

,2m  y 


(16) 


where 


a  ? 

Za 

T  - -j 

a  -  S 


and  tne  internal  eigenvectors  are: 


£■  -  r’"1  jVTu+Tf  5lra*?ta 


(17) 
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If  the  dislocation  is  located  on  the  Z  axis,  then  the 
calculation  simplifies  greatly  as  all  terms  with  m  >  2  vanish  We 
need  to  evaluate  the  functions 

l  «&>.  I  £■).  I  <&>  «  •  «  0.  *  ’  0.  r  -  L 

After  some  algeora  and  using  simplifying  relations  from  Ben-Mtnahem 
and  Singh,  1968,  r:  431),  we  find  that  the  stresses  associated  with 
a  vector  of  the  form 


u  *  ♦  f2  vtU  *  1)  3^  ♦  f.  *  1)  C. 


2-m 


(18) 


can  be  written 


Uw 


f  f 

—  *  i  — 
r  r 


7  li'-  111  4 2m 


,/f2  ,  M  1  1  (Z  *  2):  . 

-  \  —  '  —j  J 7  -  s2o 

ti.  *  r  +  i  r~— I  1  (z  *  * 

r  :r  r  r  dr  r  7  (Z  -  1 J 1  slm 


(18) 


<!±+<  rflr  ±-  II 


i  (z  *  i): 

7  rrr-rn 


'lm 


f2  f3 


iz  »  2): 

U  -  2>: 


2n 


For  the  static  eigenvectors  of  interest,  the  functions  are: 
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Jini 

f,  •  r~*T  ~  A 

^  r 

r  -i  Y  *  i  -  I 

f,  »  r  - — £ - -  z 

1  Y 

f,  .r-<1  +  2> 

4 

fx  -  -  (1  ♦  1)  r“<*  *  2> 

(20) 

1 

%  : 

f3  .  r-<  **  1) 

we  now  define  th e  coefficients  A,,,  Af,  AH  to  be  the  conplex 
conjugate  of  the  right  hand  side  of  equation  (ig)  double  dotted  with 
n  and  s,  evaluated  using  the  functions  fj  and  fj  for  the 
vectors  and  respectively  and  divided  by  the 

Green's  tensor  factors  1(1  *  1),  l  (2  i  -  1),  ( i  *  j)(2z  ♦  3) 
respectively.  The  tractions  on  the  surface  of  the  cavity  are  then: 

*■  *F  I  (f)  *  K  (21) 

*  **  I  <f  >  *  «r 

Finally,  computing  all  quantities  explicitly,  -e  get: 


dv  nl 

5U  %  ‘  “Si 


x  U  *  1 ) (2  -  v)(2£  ♦  3) 


♦  2«  U  ♦  2  -  7 )(l  *  l)2  j  5**  A.  z(£  -1)  R*'2 


Jr  n2 

Aa«  °io  ’  A| 


n  »  (2  a  ♦  i}2  -  v^  si 


(22) 


♦  Ap.  2u  (Z  -  1)  R' 
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Ali  °zm  “  Ah  “  -  *>  ** 
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APPENDIX  3 


CONVERSION  FROM  VECTOR  MULTIPOLE  COEFFICIENTS 
TO  SCALAR  MULTIPOLE  COEFFICIENTS 


In  order  to  make  synthetic  seismograms  in  layered  media,  it  is 
necessary  to  transform  from  vector  spherical  harmonics  to  vector 
cylindrical  harmonics.  The  easiest  way  to  do  this  is  to  first  form 
the  cartesian  scalar  potentials  in  sphe-ical  coordinates  and  then 
transform  these  to  cylindrical  coordinates. 


The  vector  solutions  are  written: 


"  <x*  ->  *  £  £  «*»  <“>  Jjtm  +  (“>  N-  <K.r> 


t*o  m—i 


'  M  bun  <Kar> 


1 a  '  '  J»  '  s 

(1) 


The  scalar  potentials  are  written: 

j  l»  *  u^  i  -  1,  2,  3 


X1  -  \  (v  X  uK 


X*  «  7  •  u 


(2) 


These  are  all  harmonic  functions  and  can  be  expanded  in  the  fora: 

**  ’  A  **  ■*»  (3) 

l  *0  cw- 1 


“here  .  h^2)  (K,r)  P^  (coso)  e1 


Using  the  definitions: 
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"SB 


H  *  7  x  re  *  " 
J.n  _  r  la 


(4) 


V  •  *s  7  * 

The  problem  is  to  determine  the  scalar  coefficients  in 
terms  of  the  vector  coefficients  ajm,  sZn,  vlo.  This  is 
easily  done  for  X4. 


■  7  •  u 


t,m 


’la 


*  i '  %r 

i  ,sn  a 
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i  ,n 


'rfis? 
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2a 


Comparing  with  equation  (3),  the  coefficient  is: 

A^  »  -  k  y 
tm  a  Z® 

Computation  of  x'for  i  -  1,  3  is  more  difficult. 

Z® 


(5) 


(6) 


"  T~  ?  3W  *  8tm  *£b 
im  * 

where  the  fact  the  7  •  M  «  0  has  been  used.  X*  then  are  the 
cartesian  components  of  equation  (6).  Morse  and  Feshbach  (1953) 
give  the  cartesian  components  of  the  vectors  _f^ra,  B,o  and  Clg  which 
can  be  used  to  reduce  Equation  5  to  scalar  eigenfunctions: 
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Jui  *  1)  Cl51  -  |  1  (i  ♦  <11  )(i  -  m  *  l)i  Y^"1  +  i  Y^1  J 


-jj  (i  +  b)(  t  -  a  *■  1)  y'J'1  -  Y^"1 1 


(7) 


-  k  1m  y“£ 


Using  *  '111  1*1)  and  comparing  with  equation  (3),  we  find: 

f.AL  •  8im+i  ~r  < 1  *  ra  *  UU  -  n)  ♦  a£nl_1 

Z  ~^a  Sc 

BAtm  '  3:mrl  T  +  °  *  DU  -  ns)  *  a,^  -£ 
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a  in 
Similarly, 


A'-  *  -  in  slm  T 
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Lfci 
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[tin 
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where: 


_&n 


U  ♦  a)(t  *  a  -  1}  Y®".1  -  Y^.1 

t-1  £-1 


U  -  m  *  l)(£  -0*2)  Y^|  ♦  Y^J 


*7  j 


(1  ♦  a){i  *  a  -  1)  i  Y?"}  ♦  1  yT,1 

Z--1  £-1 


-  (*  -  »  *  UU  -  a*J)  it*!-!  rf,1 
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A.2  ,-4 
o  *.m  4 


(12) 

(t  *  2) 

vitM  trrr 


••.-!,«  fei1] 

X-r  [ 


ai+i,  m 


<1  ♦  2  Hi  +  m  +  :' 

jrn - 


1-1,01  Zl  -1 


coefficients°nthat,|r<1  (I2)  ““  be  Mded  for  the  COTD,ete  "wltipole 


*  X  *  X’  1  -  1.  2.  3  (13) 


This  completes  the  transformations.  However,  it  is  sometimes 
convenient  to  nave  the  coefficients  exoressed  in  the  form: 

»  z 

x  *  £  Z)  *lm  plm  (cos»)  Jr)  cosmo 
l*o  m*o 


*  bL  f>lm  (cos®>  hi2)  <«•>  Sira# 
Using  the  relation 

pr<c°s»)  ■  i-l)a  -{f-:  «f(cos*) 


(14) 


(15) 


We  get 
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+  (-l)m  (£  -  m)i  A 

im  tr^mi'r Ai 


Bl”  %-(-1)ra|r^|tA^ 


Minster  (1976)  gives  the  transformations  for  rotations  of  the 
multipole  coefficients.  The  coefficients  derived  for  stress 
concentrations  along  the  z  axis  can  be  rotated  to  any  other 
location. 

The  displacements  may  be  recovered  from  the  potentials  using 
the  relation: 


U  (x,  «)  -  -  4-  7  X4  ♦  -^X-  7  x  x 

-  -  IT  -  - 

a  S 


Far  field  P  waves  have  the  simple  form: 


Up  -  £  I1  L  (A*  (^)  cosnip  ♦  6A  (u)  sinmp)  P  (coso) 

K  r  l  m  "  ™ 


The  scalar  coefficients  of  equation  (16)  are  now  in  the  form 
used  by  Bache  and  Karkrider  (1976)  for  computation  of  synthetic  body 
wave  seismograms  and  by  Harkrider  and  Archamoeau  (1S81)  for 
computation  of  synthetic  surface  wave  seismograms.  These  methods 
are  used  for  all  of  the  examples  in  this  report. 


SrsriMJL  tClXHCt  AHO  scrrwAAt 


APPENDIX  4 


SCALING  RELATIONS  FOR  MULTIPOLE  COEFFICIENTS 


The  displacement  field  for  tectonic  release  from  a  cavity  of 
radius  R  can  be  written 

u  («,  r;R)  .  R:  E  A^  (R>  f£  UR)  g,  («r) 

L,m 

where  A  depends  an  the  direction  of  observation  and  the  type  of 
prestress  field,  f  is  a  function  of  cavity  radius  and  g  is  a 
propagation  term.  If  ye  change  the  cavity  radius  to  R'  leaving 
elastic  constants  the  same,  then 


u  (»,  r;  R* )  .  (R* )2 


2  *, 


(R-)  f tW)  gzUr) 


If  we  look  at  the  displacement  at  frequency  u'  «  "£,  then 

u  u* .  r,  RM  .  (R')2  Z  \a  <R')  *L  (“R)  9t  (»  Jr) 

For  the  case  of  constant  prestress  fields  A  (R • )  •>  A  (R) 

2,ro  ' 

and  we  have  the  simple  relation 


In  the  far  field  g^  is  proportional  to  IUr.  Using 
1  R* 
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u  {*•,  r;ft* )  .  u  («,  r;  R) 


so  the  displacement  varies  as  the  cube  of  the  radius  at  a  frequency 
which  vari.'s  inversely  with  radius. 

Similarly,  for  the  inhomogeneous  stress  fields  considered,  the 
coefficients  have  the  form 


A2m  (R,)  *  ~S  FZa 


D 

If  • ie  vary  the  stress  field  so  that  £  -  constant,  then 

u(»‘,  r  ;R 1 )  «  ( ^  £rj  u  (u,  rj  R)  »  u  (  u,  r;  R) 

If,  however,  we  use  the  maximum  stress  (or  average  stress)  as  a 
parameter,  the  scaling  for  homogeneous  and  inhomgeneous  stress 
fields  becomes  the  same. 


For  homogeneous  stress  fields  the  displacement  is  proportional 
to  the  prestress.  For  the  stress  concentrations,  the  stress  field 
is  proportional  to  l/l3.  If  R/L  is  constant,  then  the  stress 
field  at  all  points  on  the  surface  scales  as  l“  .  So,  using  the 
stress  (at  any  point)  a  as  a  parameter,  the  scaling  relation  for  all 
of  the  stress  fields  becomes 


•’* .  a 


*a 

'  i 
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At  low  frequencies  the  spectrum  Is  approximately  flat  so  the 
displacement  (and  the  generation  of  surface  waves)  are  proportional 
to  oR  .  At  high  frequencies  the  spectrum  is  proportional  to  u 
so  the  emplacements  scale  with  «R. 


